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Guaranteed Cost Dynamic Coherent Control for Uncertain Quantum 

Systems 

Chengdi Xiang, Ian R. Petersen and Daoyi Dong 


Abstract —This paper concerns a class of uncertain linear 
quantum systems subject to quadratic perturbations in the 
system Hamiltonian. A small gain approach is used to evaluate 
the performance of the given quantum system. In order to 
get improved control performance, we propose two methods 
to design a coherent controller for the system. One is to 
formulate a static quantum controller by adding a controller 
Hamiltonian to the given system, and the other is to build 
a dynamic quantum controller which is directly coupled to 
the given system. Both controller design methods are given 
in terms of LMIs and a non-convex equality. Hence, a rank 
constrained LMI method is used as a numerical procedure. 
An illustrative example is given to demonstrate the proposed 
methods and also to make a performance comparison with 
different controller design methods. Results show that for 
the same uncertain quantum system, the dynamic quantum 
controller can offer an improvement in performance over the 
static quantum controller. 

I. Introduction 

In recent years, there has been considerable interests fo¬ 
cusing on quantum feedback control due to its applications in 
metrology, quantum optics, quantum computation, quantum 
communication and other quantum technologies [1]-[16]. 
Moreover, it has been recognized that quantum feedback con¬ 
trol plays a vital role in manipulating a quantum mechanical 
system to achieve some pre-required closed loop properties 
such as stability [1], [2], robustness [3], [4], entanglement 
[5] or other performance requirement [6], [7]. In particular, 
linear quantum optics are widely studied in control areas and 
a quantum optical system can often be described by a set of 
linear quantum stochastic differential equations (QSDEs) [2]. 

In the conventional picture of quantum feedback control, 
digital or analog electronic devices are often implemented 
as controllers [8]. However, this tends to destroy quantum 
coherence and involves the destruction of quantum infor¬ 
mation in the process of making measurements. Hence, 
recent research has been focused on using a fully quantum 
system as a controller, which is referred to as a coherent 
controller, e.g., [6], [9]. Compared to measurement feedback 
control, the advantages of coherent control are to preserve 
quantum coherence, to achieve improved performance and to 
obtain a high speed processing bandwidth. However, when 
we consider coherent feedback control using the QSDE 
description, the issue of physical realizability of controllers 
arises [2], [9]. That is, in QSDEs framework, the state 
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space matrices defining the coherent controller are required 
to satisfy certain conditions in order that the controller 
represents a physically meaningful quantum system and we 
call these kind of conditions physical realizability conditions. 
However, in this paper, we use an (S', L, H) framework to 
define quantum systems, where H is a Hamiltonian operator, 
L is a vector of coupling operators and S is a scattering 
matrix [10], [11]. The quantity L describes the interface 
between the system and the field, and the operator H defines 
the self-energy of the system. Since the parameters (S, L, H) 
already describe a physically realizable system, we do not 
need to be concerned about physical realizability conditions. 

Quantum feedback controllers have been designed with 
a number of different techniques. Eor example, an H°° 
synthesis approach [2] and a LQG method [9] have been used 
to design a quantum controller for a class of linear quantum 
stochastic systems; [14] used a transfer function method 
to analyze the robustness of feedback quantum systems. 
Nevertheless, few papers have considered quantum controller 
design based on an {S,L,H) description. Based on the pa¬ 
rameters {S,L,H), we are going to design a guaranteed cost 
coherent controller not only to robustly stabilize the uncertain 
quantum system, but also to guarantee a specific level of 
performance for any admissible value of the uncertainties. 

In the previous papers on quantum controller design [2], 
[9], the coupling between the plant and the controller is 
via a field coupling which we call indirect coupling. In the 
controller design parts of this paper, we use two different 
methods. One is to add controller Hamiltonian and the other 
is to construct a directly coupled quantum controller for 
the given system. Here, direct coupling refers to that two 
independent quantum systems may interact by exchanging 
energy [15], [16] and this energy exchange is often described 
by an interaction Hamiltonian. 

In this paper, the guaranteed cost coherent controller 
design is given in terms of LMI and nonlinear equality 
conditions. We use a nonlinear change of variables to convert 
the problem into a rank constrained LMI problem which can 
be solved using an alternating projections algorithm [18]. 

This paper is organized as follows. In Section |I^ we 
define the nominal quantum system under consideration as a 
linear system using parameters (S', L, H). Then an uncertain 
perturbation of the Hamiltonian is introduced in terms of a 
commutator decomposition and sector bound conditions in 
Section imi In Section |IV] the cost function for uncertain 
linear quantum systems subject to quadratic perturbation of 
the Hamiltonian is defined and a small gain type performance 
analysis result is presented. In Section we introduce a 


controller Hamiltonian and present a theorem to show the 
construction of this guaranteed cost quantum controller. In 
Section |Vl] a dynamic controller system is directly coupled 
to the uncertain quantum system. The corresponding con¬ 
troller design and numerical procedures are presented. An 
illustrative example is presented to demonstrate the coherent 
controller methods in Section VII We also make a per¬ 
formance comparison between the static coherent controller 
and the dynamic coherent controller. Some conclusions are 


presented in Section VIII 


II. System Description 


The open quantum system under consideration is an uncer¬ 
tain linear quantum system defined by parameters {S,L,H), 
where H refers to the system Hamiltonian and can be 
decomposed as Tf = Hp + Hun- Here, Hp denotes a 
known nominal Hamiltonian and Hun denotes a perturbation 
Hamiltonian contained in a specified set of Hamiltonians W 
[3]. We assume that Hp is in the form of 


Hp 


^ p'^ ]Mp 


q 

p 


( 1 ) 


To proceed, we also define the corresponding generator 
operator as follows 

0(X) = -i[X,7J]-f/:(X), (7) 

where £(X) = [X, L] + \ X]L. Here, the notation ^ 

stands for the adjoint transpose of a vector of operators and 
[X, H] = 'KH — H'K describes the commutator between two 
operators. The following lemma will be used in the main 
results presented in this paper. 

Lemma 1: [4] Consider an open quantum system defined 
by (S', L, H) and suppose there exist non-negative self- 
adjoint operators V and W on the underlying Hilbert space 
such that 

g{y) + w<\ ( 8 ) 

where A is a real number. Then for any plant state, we have 

1 r 

limsup— / {W{t))dt<X. (9) 

T— )-oo T Jo 

Here W(t) denotes the Heisenberg evolution of the operator 
W and (•) denotes quantum expectation; e.g., see [4] and 
[ 11 ]. 


where Mp is a real and symmetric Hermitian matrix with 
dimension 2np x 2np. Here, g is a vector of position 
operators and p is a vector of momentum operators. The 
commutation relations between position and momentum op¬ 
erators are described as follows 
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The coupling operator L is of the form 


L=[ N2 



(3) 


where Ni G and N 2 G We also have 


L 

L* 


Np 


q 


■ N, 

N2 

q 

P . 



N* _ 

. P . 


(4) 


We consider self-adjoint “Lyapunov” operators V in the 
following form 


V = [ g^ P^ ]X 


q 

p 


(5) 


HI. Perturbation oe the Hamiltonian 

In this section, we introduce a perturbation for the quan¬ 
tum system under consideration. First, we define the per¬ 
turbation of the Hamiltonian in terms of a commutator 
decomposition. Then, we introduce the formulation of the 
quadratic perturbation in the system Hamiltonian. 


A. Commutator Decomposition 

For the set of non-negative self-adjoint operators V and 
given real parameters 7 > 0, i5 > 0, a particular set of 
perturbation Hamiltonians Wi is defined in terms of the 
commutator decomposition 

[V,Hun] = [V,z^]w-w'^[z,V] (10) 


for V GV, where w and 2 : are given real vectors of operators. 
Wi is then defined in terms of sector bound condition; 

vFw<^z^z + 5. ( 11 ) 

T 


We define 


Wi = 


Hun : 3 such that (101 and (11) 
are satisfied V VG V 


( 12 ) 


Lemma 2: Consider an open quantum system {S,L,H) 
where H = Hp + Hun and Hun G Wi, and the set of non¬ 
negative self-adjoint operators V. If there exists a V G V 
and a real constant A > 0 such that 


where X G K 2 Ti_px 2 np ^ symmetric positive definite 
matrix. 

Therefore, we define a set of non-negative self-adjoint 
operators V as follows 

V of the form (j^ where X > 0 is a 1 
symmetric positive definite matrix j ' 


-i[V, Hp] +C{V) + [y, 2 ^] [ 2 , V] + ^z'^z + W <X, (13) 
then ^ 

limsup;(|(; f {W{t))dt < X +6,yt>0. (14) 

T-j-oo T Jq 




































(16) 


Proof: As we know V G V and Hun G Wi, 
g{V) = -i[V, Hp] + C{V) - i\V, z^]w + iw'^[z, y], (15) 
Since V is symmetric = [z, y]. Therefore, 

0 < ([y z'^] - iw^)([y, 

= [y, z^][z, y] + r[y, z'^]w — ivF[z, V] + w'^w. 

Substituting into GD and using the sector bound 

condition ( |TT| , the following inequality is obtained: 

G{V) < -r[y, Hp] +C{V) + [y, z^] [z, y] + ^z'^z + 5. 

7 

(17) 

It follows from that Q(V) + W < A + <5. Consequently, 
the result follows from Lemma [T] □ 

B. Quadratic Hamiltonian Perturbation 

A set of quadratic perturbation uncertainties is dehned in 
the following form 


1 


Hun = ^[ Ci Q 


C, 

Cp 


1 ^ 

where = A. We also have the relationship 

= E 


C, 

Cp 


Hence, we write 


1 


Hun = ^[ P^ ] E'^AE 


The matrix A is subject to the norm bound 

l|A||<^ 


yy / Hun of the form ( [T8] l such that 


condition ( |2T] i is satished 
Lemma 3: For any set of self-adjoint operators V, 
m c Wi 

Proof: Given any Hun G W 2 , let 


1 

^=2 


A?; 


A 


12 


c, 

A 22 J L Cp. 


= 2 ^^ 


and 


Hence, 


Also, 

Tr 

W [z 

Hence, 


1 / CjAnC,V^ + CjAi2Cpy \ 

2 \ +Cp Af2Cqy + cjA22Cpy / 

-I / J.'T T ^ A M . ^'T-r T- A 


CjyAnC, + Cr^Ai2Cp 

+Cj ^ Af 2 Cg + A 22 CP 


( 28 ) 


[y,z^]rr;-'u;'^[z,y] 

_ 1 / yCjAnC, + f"CjAi2Cp 
2 V +i^cjAf2C9 + ycjA22Cp 
_ 1 / CjAnC,f" + CjAi2Cpf" 

2 V ECpAf2(qV + (pA22(pV 
= VHun-HunV=[V,Hun] 


(29) 


Also, 

1 


- c c 

4 L ^9 


AA 


c, 

Cp 


CJ] 


C, 

Cp 


(18) Therefore, we have >y 2 C Wr- 


(30) 

□ 


(19) 


( 20 ) 


( 21 ) 


IV. Performance Analysis 

In this section, we evaluate the performance of the given 
uncertain linear quantum system. First, we need to define the 
associated cost function for a quantum system as 

fT 

)dt (31) 


J = limsup;^ [ ([ p'^ ]R 

T—s-oo J Jo 


where R > 0. We denote that 


W = [ q'^ p'^ 


R 


(32) 


where ||.|| refers to the matrix induced norm and we have 
the following definition 


( 22 ) 

(23) 

(24) 


In order to introduce our results on performance analysis, 
we require the following algebraic identities. 

Lemma 4: Consider V G V, Hp is of the form Q and L 
is of the form Q. Then we have 

n T 


[V,Hp] = 


and 


/:(y) = -2 


{XEMp - MpEX) 


1 1 

= E 

q 

(25) r 


. P . 

where J = 


■ Tr(AEA|, 


I 0 
0 -I 


{nIjNp^X + X^nIjNp) 

NpE), 

. Moreover, 


I 0 
0 0 


(33) 


q 

p 

(34) 


Hun = W^z = -[ ] E'^AE 

Then, for any V G V, 

_1 ( VC^A,,Cq + VC^A,2Cp 
\^-2[+VCjAUq + VC^A22Cp 
_\( cJf^AiiC, + cJ^yAi2Cp 
2 V +Cjl^Af2C9 + CJ1^A22 Cp 


(26) 


A 


= 2SA 


(35) 


Proof: The proof of these identities follows via straightfor¬ 
ward but tedious calculations using Q. □ 


(27) 


Lemma 5: For V G V and z dehned in ( [T^ , 
[z, y] = 2E^X 


( 36 ) 




















































[F,z^][z,F] =4 


1 T 




T 

Z Z = 


1 T 


E^E 


(37) 


(38) 

□ 


Proof: The result follows from Lemma 

We now in a position to present our main result in this 
section. 

Theorem 1: Consider an uncertain quantum system 
(S', L, E[), where El = Elp + Hp is in the form 

of (|T]i, L is of the form 0 and G W 2 . If Gl = 

—i'SMp — f EiVp JiVp is Hurwitz, and 


A^X + XA+ ^ + R 2XJ:'<E^ 
2ET,X -I 

has a symmetric solution X > 0, then 

1 

J = limsup— / {W{t))dt 

T—>00 7 Jq 

= limsup ^ f {[ q'^ p'^ ] R 

T-J-OO r Jq 


< 0 (39) 


)dt < A + 5 


where 


A = TifXENl 


I 0 
0 0 


NpE). 


(40) 


(41) 


Proof: The proof is similar to that in Theorem 1 of [6], and 
a detailed proof is omitted. □ 

V. Static Coherent Controller 

In this section, we aim to design a coherent guaranteed 
cost controller for the given uncertain quantum system by 
adding a controller Hamiltonian EIki- As there are no 
additional dynamic variables defined by this Hamiltonian, 
we refer this kind of quantum controller as a static coherent 
controller. 

The controller Hamiltonian EIki is assumed to be in the 
form 


Hki = 1 


qT pT 1 


(42) 


where F G K 2 npx 2 np ^ 


X2np : 


IS a symmetric 


matrix. An associated cost function J is defined in the 
following form 


J = lim sup — 

T—^OO 


1 T 


(R+pF^KFE^KF) 


where p G (0,c») is a weighting factor. We have 


W = 


-I T 


(R + pF'^KFF'^KF) 


(43) 


(44) 


A. Controller Design 

Theorem 2: Consider an uncertain quantum system 
{S,L,H), where H = Hp + H^n + Hki^ Hp is in the 
form of Q, L is of the form (|^ and Hun G W 2 , and the 
controller Hamiltonian Hki is in the form of ( |42l l. If there 
exists symmetric matrices K, X > 0 and Y — KFO'^X 
such that 

AXe^E^ F^KF 

-I 0 <0 (45) 

0 -IIP 


B 

AEOX 

F'^KF 


where B = A^X + XA + 2F'^Y + 2Y'^F + E'^E/-f‘^ + R 
and A = —lYMp — ^YNpJNp, then the associated cost 
function satisfies the bound 


1 

limsup — / {W{t))dt < ^ + 5, 

T—loo 7 Jq 

where W (t) is in the form of (®| and 


A = TriXENj, 


I 0 
0 0 


Proof: Suppose the inequality (45i is satisfied. Using the 


NpY). 


(46) 


(47) 


Schur complement [19], we have 

B + IQXe'^E'^EOX F'^KF 
F^KF -Ijp 


< 0 . 


(48) 


Also, it follows from the Schur complement that (481 is 
equivalent to the following inequality 

A^X + XA + 2F^Y + 2Y^F + E^E/j^ 


+ leXe'^E^EOX + R + pF^ KFF^ KF < 0. 


(49) 


Substituting Y = KFO^X into the inequality (49 1 , we 
obtain that 

{A - iYF'^KF fX + X{A - iYF^KF) 

+ AXY'^E^EYX + E'^Ej-i^ ere pF'^KFF^KF < 0. 

(50) 


It follows straightforwardly from (50i that A — iYF'^KF is 
Hurwitz. Based on Lemma and Lemma we have 

- ^[U, Hp E Hki] + C{V) E [U, z^Uz, U] + ^ 

7 


- Tr(XSNl 


/ {A-iYF'^KFfX 
EX{A-iYF^KF) 
EAXY^E'^EYX E E'^E/-f^ 
\ ere pF^KFF'^KF 

' ol«pS). 


}dt, 


(51) 

It follows from Lemma and Lemma that the following 
conditions are satisfied; 

limsup f {W{t))dt < X E S, (52) 

T-100 T Jq 


where A = Tr(XSAf], 


10 

I 0 

0 0 


NpY). 


□ 














































Remark 1: In order to construct the guaranteed cost dy¬ 
namic controller in Theorem]^ we should not only satisfy the 
linear matrix inequality ( |45] ), but also guarantee the equality 
Y = KFO"^ X. We should notice that the above equality is 
non-convex constraint. Hence, we will convert this problem 
into a rank constrained LMI problem. The numerical method 
is similar as that in section [Vij-B. Thus a detailed procedure 
is omitted in this section. 


VI. Dynamic Coherent Controller 


In this section, we design a dynamic coherent controller 
that is realized by being directly coupled to the original 
uncertain quantum system. The directly coupled controller 
and the quantum plant may interact with each other by ex¬ 
changing energy via an interaction Hamiltonian. We should 
also note that the associated operators of the nominal system 
commute with all the associated operators of the dynamic 
controller. In the following part, we will introduce the 
design of dynamic quantum controllers and then formulate 
the augmented quantum system consisting of the quantum 
plant and the directly coupled quantum controller using an 
{S,L,H) parameterization. 

We consider a linear dynamic quantum controller Hamil¬ 
tonian in the form of 


Hk2 = ^ 



T 

P -l 

qx 

K 22 

qx 

Pk 

Pk 


(53) 


where K 22 S K 2 njfx 2 nx follows, we take the 

interaction Hamiltonian to be 

T 


_1 
Hint —77 


Pk 


KI 2 F 


F^Ki2 


QK 

Pk 


(54) 


where K 12 G '^'^■npx2nK^ Also, we consider a static coherent 
controller in the form of ( |42l l. 

It follows from ( |42| i, ( |53| ) and ( [54l i that the total controller 
Hamiltonian is in the form of 


Hk = 

_ 1 
“ 2 


Here, we denote K = 


Hki + Hk2 


P _ 

T 


- -| 

q 


q 

p 


F'^KiiF F'^Ki2 


p 

qx 


Kl^F K 22 


qx 



. Pk 


(55) 


Kii 

KTo 


Ki2 

K 22 


_ 

■ F 

0 ' 

, and F = 

0 
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Consequently, the total Hamiltonian for the augmented quan¬ 
tum linear system (excluding the perturbation Hamiltonian) 
that comprises the nominal quantum plant and the quantum 
controller is 

Hp+k = Hp+Hk. (56) 


We know that the coupling operator Lp and 


Lp 


L 


# 


for 


the plant are defined in Q and Q, respectively. The coupling 


operator for the controller is assumed to be in the form 


Lk = [ Hxi Nk2 ] 


Qk 

Pk 


( 57 ) 


where Nki G and Nk 2 G 

Also, we have 

Hki Hk2 


Lk 

r# 

= Nk 

qx 

_ 

. Lji 


. Pk 



^^Kl 


N 


# 


K2 J 



qx 


. Pk . 


(58) 

Hence, the overall coupling operator for the directly coupled 
system comprising the plant and the controller is of the form 


Lp 


■ iVpi 

Np2 

0 

0 

Lk 


0 

0 

Nki 

Nk2 


L = 


Also, we write 


L 

L* 


q 

p 

qx 

Pk 


(59) 


Lp 
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Lk 

T* 

Lp 

= N 

p 

qx 


r# 

^K J 


. Pk 


Wpi 

Np2 

0 

0 

0 

0 

Nki 

Nk2 

JVpi 

AT# 

-''P2 

0 

0 

0 

0 

^^Kl 

^^K2 
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qx 


. Pk 


(60) 


The “Lyapunov” operator V for the augmented system is 
considered to be a non-negative self-adjoint operator of the 
form 


V = 


where X = 


P ^ 

T 


p -l 

q 


q 

p 


Nil Ni2 


p 

qx 


Ni 2 N 22 


qx 

. Pk 


. Pk 


( 61 ) 




^12 

^22 


(np+np:)x 2(np+nj^) 


IS a 


symmetric positive definite matrix. 

Hence, a set of non-negative self-adjoint operators P can 
be defined as 


V = 


V of the form (61) where X > 0 is a 


symmetric positive definite matrix 


(62) 


Since all variables of the plant are assumed to commute 
with all variables of the controller, we have the following 
commutation relation 
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q 

T - 
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p 



qx 
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qx 


_ 

. Pk 


. Pk 

_ 


= 2ie = s 


(63) 


where 0 = 


0 


In order to proceed to the 

following section, we introduce the permutation matrix 
Pn+m, where the symbol Pn+m refers to a 2{n + m) x 
2(n -f m) matrix. The permutation matrix Pn+m is de¬ 
fined in such way that if we consider a column vector 





































































a = [ai 02 ... 
[ai 02 ... a„ 
^2n-\-n 


^n+m ■■■ ^2(n+m)] 5 then — 

t^n+m tl2n+m+l 


notation we need to introduce is = 
we have 


0 -In 


. Hence, 


Pr 


mp-\-mK 


N = 


N = P, 


Pi 


mp+niK 


mp+iriK'^'^P+^K ^rnp+niK 


Np 0 
0 Nk 

Np 0 

0 Nk 

Prr 


(64) 


'Trip 
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Jn 


In order to present the main result, we require some 
algebraic identities. 


Lemma 6: Suppose V GV, Hp+k is in the form of (561 
and L is in the form of ( |60| ). Then, we have that 

-i[V,Hp+K]+CiV) 


(65) 


q 

T 

/ T' \ 

q 

p 

( {A-iEF KFfX \ 

P 

qK 

\ +X{A-iEF'^KF) j 

qK 

. Pk 


. Pk 


+ Tr(X57Vi 


I 0 
0 0 


NE), 


^n+1 ^n+2 

fl 2 (ra+m)]^- Recall the property of an m x m permutation 
matrix that it is a full-rank real matrix whose columns 
comprise standard basis vector for M™, that is, vectors in 
K™ contains precisely a single element with value 1 and 
all the remaining elements are 0. A permutation matrix P 
also has the unitary property PP^ = P^P = I. Another 

'In O' 


Following (64 1 , we get 

Mp + F^KiiF F^N 


A = — 

1 , 

~ 2^ 


12 


K^F 
■ nI 0 
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22 


0 
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mp 

0 J, 
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JniK 

-iUMp - IHN^JNp 0 


Np 
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0 

Nk 
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- iEF^KF 
=A - iEF^KF. 


-Ij:nI,jnk 


(69) 

Therefore, the result follows from ( |67l ) and ( |69l l. □ 

We know that only the nominal quantum system is sub¬ 
ject to uncertain quantum perturbation. Hence, we rewrite 
the quadratic perturbation uncertainty ( p0| ) and in the 
following form 
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QK 

Pk 


1 ^ 

q 

[ E 0 A [ A 0 ] 

p 

qK 

. Pk 


(70) 


z = 


c 

Cp 


« =[ E 0 ][ q'^ qj. ] . (71) 


In order to better present the following theorem, we denote 
E=[ E 0 ]. 

Now we are in the position to define an associated cost 
function as 


1 r°° 

J = limsup — / {W)dt 
T—>00 I Jo 


and 


(72) 


where A = 


-iY,Mp - iSWpWp 




'‘K 


Proof: According to Lemma we know that 
-i[V,Hp+k]+ C{V) 


( 66 ) 
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f {-iEM-\ENfjN fX \ 
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+X{-iEM - \ENUN) ) 
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■ TriXEN'^ 
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TVS), 


where M = 


Mp + E'^KiiF E^Ki2 
K 22 


KI 2 F 


(67) 


To shorten the 


writing, we denote A = —iEM — ^SA^ JA and we have 


A = — iS 


Mp + F^KiiF F^Ki2 

KI 2 F K 22 


( 68 ) 


- -EN'fJmp+rUKN 
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qK 
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qK 

Pk 


pe'^kef'^kf 
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qK 

Pk 


1 ^ 
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(A -f pF'^KFF^KF) 
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qK 


. Pk . 


R 0 
0 0 


(73) 


and p G (0, 00 ) is a weighting factor. 


A. Controller Design 

In this subsection, we present the main result on dynamic 
coherent controller design. 

Theorem 3: Consider an uncertain quantum system 
{Sp,Lp,Hppnn), where iTp+„„ = Hp + H^n, Hp is in 
the form of (0, Lp is of the form ([^ and Hnu G W 2 . 
It is directly coupled with a dynamic quantum controller 
{Sk, Lk, Hk)- Then the closed loop quantum system is 
defined by {S,L,II) where H = Hp + Hun + Hk, Hk 
is in the form of (l55ll and L is in the form of (|59ll. If there 



























































exists a symmetric matrix K, X > 0 and Y = KFQ^X 
such that 


B 

AEQX 

—T — 

F KF 

-^T , 


4X0'^^;^ F^KF 


-I 
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0 

-I/p 


< 0 


(74) 


where B = A" X + XA + 2F" Y + 2Y'^F + E E/^"^ + R 
and A is in the form of ( 661 , then the associated cost function 
satisfies the bound 

limsup;^ f {W{t))dt < X + S, (75) 

T->-oo T Jq 


where W (f) is in the form of (73 i and 
A = Tr{XEN^ 


I 0 
0 0 


iVE 


(76) 


Proof: Suppose the conditions of the theorem are satisfied. 
Using the Schur complement, 0 is equivalent to 


B + IQXQ'^E^EQX f'^KF 


—T — 

F KF 


-I/p 


< 0 . 


Applying the Schur complement to ( |77] l, we obtain 

A^ X + XA + 2f'^Y + 2Y^F + e'^E/j^ 

+ leXe^E^EOX + R + pP'^KFE^KF < 0. 


(77) 


(78) 


Substituting Y = KEQ^X into (78i, we have 

(A - iEF'^KF fX + X{A- iEF^KF) 

+ AXE^lfEEX + lfE/-f‘^ FR + pF^KWp'^KF < 0 




(79) 

T ■ 


It follows straightforwardly from ( |79| that A — i^F KF is 
Hurwitz. Following Lemma (|70| and (|7T]i, we have 


[V,z ][z,U] + -z ^ 


7 

P 

Qk 

Pk 


ixe^e'^eex 

+lfE/-f^ 


q 

p 

Qk 

Pk 


(80) 


Based on (80i and Lemma the following equation is 
achieved; 


- i[V, Hp+k] + C{V) + [U, z^] [z, Vj + ^z^z + W 

7 

/ (A-iEF^KF)'^X \ 

+X{A - iEP’^KF) 
FAXE'^E^EEX 
+E^E/-f^ + R 

—T ^.-F=rF^T 


q 

T 

p 


qK 


. Pi< 
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qK 


. Pk 


+ Tr(A:SAr1' 


I 0 
0 0 


TVS). 


( 81 ) 


According to Lemma and Lemma we obtain 


lim sup — 

T —¥oo ^ 


where A = Tr(XSAf'^ 


0 

/ 0 
0 0 


{W{t))dt < A + (5, 
TVS). 


( 82 ) 


□ 


B. Numerical Procedure 

Theorem 13 can be formulated as the minimization of the 


cost bound A subject to the LMI (74i and the equality 
Y — KFQ^X constraint. Due to a fact that the equality 
constraint is not linear, we cannot use a traditional LMI 
technique. Therefore, this problem is transformed into a 
rank constrained LMI problem [18]. The way to deal with 
this problem is to linearize the equality Y = KFQ^X 
by introducing appropriate matrix lifting variables and the 
associated equality constraints, and transforming them into 
an LMI with a rank constraint. 

We set Z^i = X,Zx 2 = K and introduce appropriate 
matrix lifting variables Z^i = KFQP", Zy 2 = KFQ^X. 
Then we define Z as a symmetric matrix with dimension 
5(2np + 2np-) x 5{2np + 2np-) and Z = VV^, where 
V=[I zl, zl^ zji 

We require the symmetric matrix Z to satisfy the following 
conditions: 


Z > 0; Zop — 1 = 0', 

Zxi — Z'^i= 0 ; Zx2 — z'^2 = 0 ; 

Zvi ~ Zx2FQ"^ = 0; Z„2 — ZyiZxi = 0 


(83) 


and a rank constraint 


rank(Z) < 2np + 2np;- (84) 


Therefore, the way to realize Theorem is to solve ( fT^ I 
subject to condition ( [S^ l and ( (84) > using a rank constrained 
LMI. To solve the rank constrained LMI problem above, 
we work in Matlab using an algorithm [17] and also use 
the Yalmip optimization prototyping environment [21] and 
SeDuMi solver [22]. 

In order to minimize the cost bound A, we also need to 
introduce another inequality 


TviXEN"! 


I 0 
0 0 


TVS) < 13. 


(85) 


Hence, to minimize the cost bound, we just need to lower 
the pre-specified bound /3 using the bisection method; e.g., 
[ 20 ]. 


VII. Illustrative Example 


In order to illustrate the controller design methods demon¬ 
strated in this paper, we are now present an example. We 
consider an uncertain quantum system defined by a triple 
{S,L,II) in the following form 


Sp = I,Lp = ^ [ 
Kp-\-un ~ 2 ^ ^ 



'0.7 1 ■ 
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1 1.8 


. p . 



































where Hp+un = Hp+Hun- The Hamiltonian of the nominal 
system is 




T T 

q p 
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. P . 


and the Hamiltonian of the perturbation system is 
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. P . 


This leads to the corresponding parameters in Theorem 
and Theorem |3] as follows: 


Mp = 


Nk = - 

2 
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VO-di 
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-\/Ki 
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1 

,A = 


By running the program, we find that the nominal plant 
is mean square stable for the range k G (8.2, c») without 
any controller. Nevertheless, a static coherent controller can 
stabilize the system for the range of k € (7.1, oo). Mean¬ 
while, a dynamic coherent controller can guarantee that the 
system is stable for the range of k G (7, oo). Therefore, we 
can conclude that quantum controllers can stabilize uncertain 
quantum systems for a larger range of k than the one without 
a controller. 

With regard to the stability and the performance of the 
quantum systems, the importance of the static controller to 
the given system has already been addressed in [6]. Here, 
we will focus our attention on a performance comparison 
between the given system with a static quantum controller 
and that with a dynamic quantum controller. As displayed 
in Figure 1, compared to the static controller, the dynamic 
quantum controller can guarantee the uncertain quantum 
system stable for a larger range of k values with a lower 
cost. Hence, it leads to a closed loop system having better 
performance. 


VIII. Conclusions 

This paper has evaluated the performance of uncertain 
linear quantum systems subject to quadratic perturbations 
in the system Hamiltonian. We proposed methods for a 
static coherent quantum controller design and a dynamic 
coherent quantum controller design. We converted the con¬ 
troller design problem into a rank constrained LMI problem 
and solved the problem based on an alternating projections 
algorithm. An example was presented to illustrate these 
two controller design methods. By comparison, we showed 
that a dynamic controller can have an improved control 
performance over a static one. 
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